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Recent studies on the low dimensional systems have renewed great interests in the 
Gutzwiller-Jastrow wavef unctions. The wavefunctions are useful in the aspects that they 
may serve as good variational wavefunctions or they may be exact solutions of the Hamilto- 
nians p|^T9|. Most recently, the intense research in the field of integrable systems has shown 



the wavefunctions to be exact solutions of some quantum many particle systems. These 
integrable systems are characterized by that the full spectrum may even be written in terms 
of more generalized Jastrow wavefunctions, as in the cases of Fermi or Bose gases, 
Haldane-Shastry spin chain, and the ID supersymmetric t-J model of hopping and 
exchange. 

Hubbard model has been of great interest since the discovery of the high Tc supercon- 
ductivity. About two years ago, Gebhard and Ruckenstein introduced the one dimensional 



SU{2) Hubbard model of 1/r hopping and on-site energy U [|rT|. The model is completely 
integrable for arbitrary on-site energy. In the strong interaction limit U = oo, it has been 
discovered recently that a set of Gutzwiller-Jastrow wavefunctions is exact eigen-functions 
of the Hamiltonian the system exhibits the spin-charge separation in the full Hilbert 
space [0,M. 



In this work, we introduce a new integrable model, the one dimensional 1/r multi- 
component Hubbard model. In the following we only discuss the strong interaction limit 
U = oo. Generalizing our previous work, we show that a set of SU{N) Gutzwiller-Jastrow 
wavefunctions is eigenstates of the system. The full excitation spectrum and the thermody- 
namics are also given explicitly in this strong interaction case. Spin and charge are decoupled 
in the full Hilbert space, the system is a special example of SU{N) Luttinger liquids. At 
the end of the work, we also discuss speculations of further investigation of the system of 
finite on-site energy U. 

The Hamiltonian for the one-dimensional Hubbard model is given by 

H = Y.Y1 tiMaCja + niaUi^,, (1) 

where cl^. and Qo- are creation and annihilation operators at site i with spin component a. 



The sum over a runs from 1 to A^, where is the number of flavors of the fermions. We take 
the hopping matrix tij = it{—\)'^^^'-'^/d{i~j) where d{n) = (L/n) sin(n7r/L), and U is the on- 
site energy. Here, because of the special form of the hopping matrix, for the wavefunctions 
of the system, we assume periodic boundary conditions for odd L, or anti-periodic boundary 
condition for even L. 

In the strong interaction limit U = oo, each site can be occupied at most by one particle. 
In this limit, we work in the Hilbert space of no double occupancy and no multi-occupancy. 
The Hamiltonian can be written in terms of the Hubbard operators: 

H= E^^.^f^f- (2) 

a=l,2,---,N iytj 

Let us denote the number of holes by Q, that of the fermions of the first fiavor by Mi, 
that of the second fiavor by M2, ■ ■ ■, that of the A^— th fiavor by Mjy. Following notations 
used in previous literatures, states in the Hilbert space can be represented by spin and hole 
excitations from the state full of fermions with the A^— th fiavor \P), 

!$>= E H{^th{yj}mblllh]\p>, (3) 

Here = c\^CiN annihilates one N — th fiavored fermion at site i and creats one a — th 
fiavored fermion at site i for a = 1, 2, ■ ■ ■ , (A^ — 1), while /i] = CjN creats a hole at site j. 
Here or in the following we always implicitly assume that we work in the space of no double 
occupancy and no multi-occupancy in the discussion of the strong interaction limit. The 
amplitude {yj}) is symmetric when exchanging the fermions at positions and x", 

and antisymmetric in the hole positions {Vi}- 

Let us consider the following generalized SU (N) Gutzwiller-Jastrow wavefunctions cor- 
responding to uniform motion and magnetization [P], 

27ii 

'^Gi{xt},{ym}) = exp— [^J«^X^ + A^t/i] X $0 

^ a i i 

*o = n d'i^t - ■ n di^t - ■ n ^ixt - ■ n divm - vn) (4) 

a;i<j a</3;i,j a,i,m, m<n 

where a, (3 = 1,2, ■ ■ ■ , N — 1. The quantum numbers Ja and Jh govern the momenta of the 
fermions and the holes. They can be integers or half integers such that the wavefunctions 
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are periodic (or anti-periodic) for odd L ( for even L) under the translations xf ^ xf + L, 
and i/m — Vm + L. For the wavefunctions to be eigen-states of the Hamiltonain, the quantum 
numbers must be choosen from some restricted regions, which are to be specified below. 

To demonstrate that the wavefunctions are eigen-states of the Hamiltonian, we have to 
consider the effect of the hopping operator very carefully. The hopping operator can be 
broken into N parts, each corresponding to the hopping of fermions of different flavors. Let 
us first consider the hopping operator of the N—th fiavor, T{N) — UjC^mCjN- When it 
operates on the wavefunctions, the hopping of the fermions of N—ih. fiavor is equivalent to 
the hopping of holes, 

n=i n m^n {a,i) 

where F„^(n) = cos ^ sin ^ cot 6„^, = cos ^ sin ^ cot 6n(m); Qnm = T^iVn - 

ym)/L, ©n(m) — '^{Un — ^?)/L. The sum can be carried out after expanding the products 
and classifying terms by the number of particles involved. In the end, only the zero particle 
term and two particle terms are left. Many particle terms vanish, yielding the following 
result: 



-QJh + {2T^t/L)i J2 E cot e^iai)- (6) 



n a,', 



This result is valid under the condition \Jh\ < L/2 - [Q + (Mi M2 H h Mjv_i)]/2. 

To consider the effects of other parts of the hopping operators, we can not use the 
wavefunctions directly, since the hopping will involve the fermions and holes simultaneously 
when they operate on the wavefunctions. We can generalize the idea of spin-rotated version 
developed in the recent work of the t-J model to this SU{N) case. For example, to 
deal with the hopping operator T{N — 1) = Y^i^^j cJ(Af-^i)'^i(A^-i)' write the Gutzwiller- 

Jastrow wavefunctions in terms of the hole positions and the positions of the fermions of 
fiavors excluding the {N — l)-th fiavor. In terms of these coordinates, the wavefunctions can 
be found to be still in a similar product form, and thus effect of the operator T{N — 1) can 
be calculated in the same way as for T{N). 



For the other hopping operators T{N - 2), T{N - 3), • • • , T(2), 7(1), similar procedures 
can be carried on. After adding all the effects of the hopping operators together, the two 
particle terms vanish since positions of all the fermions and the positions of the holes span the 
entire lattice. Thus the Gutzwiller-Jastrow wavefunctions are found to be exact eigen-states 
of the Hamiltonian, with eigen-energies given by 

E{Jh- Ji, J,,..., J;v-i) = -i27rt/ L)Q[J, + J« + + . . . + jf "D], (7) 

where we have jf^ = A - Jn-i + L/2, jj^^ = A - Jn-2 + L/2, ■ ■ •, J^"^"'^ = - + L/2. 
The many particle terms vanish, and thus our result holds, under the conditions 

141 <M^/2, 
|ji')|<Mjv-i/2, 
|jf I <MAr-2/2, 



< Ml/2. (8) 

Here the ground state energy is given by Eq = — (27r|t|/L)[L/2 — Q/2]Q. 

For this multi-component system, the spectrum can also be written in terms of more 
generalized Jastrow functions. Here, we just write down the spectrum without getting into 
the detailed algebra, 

Q 

E = -{27Tt/L)J2Ki + {7TtQ/L){L + 1), (9) 

i=l 

where Ki takes values from the region (1,2, Here the result indicates that the 

spectrum is invariant when changing the sign of t. Each energy level is determined by 
a charge configuration such as (101010) for Q — 3,L = 6, where the ones represent the 
values occupied by the charge momenta K^. In this system, the spin and charge degrees are 
decoupled from each other in the entire Hilbert space. On these physical grounds, we see 
that for each charge configuration the degeneracy of the corresponding energy level is given 
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by the number of the ways to distribute the free spins among the L — Q empty values. With 
this result, we find the free energy per lattice site given by 

F(T,/x)/L = -/x-^ rdgln[iV + e^(«*-^)], (10) 

where /i is the chemical potential of the fermions. This free energy has also been found to 
correctly reproduce the first three terms in the high temperature perturbation expansion. 

In summary, we have solved the multi-component Hubbard model in the strong inter- 
action limit. In this limit, the spin degrees of freedom decouple from the charge degrees of 
freedom in the entire excitation spectrum, and the system is a special example of the SU (N) 
Luttinger liquids in the sense of Haldane. We have shown that the SU (N) Gutzwiller- 
Jastrow wavefunctions are eigen-states of the Hamiltonian. 

In the end, we notice that in the half filling and large U limit, our model reduces to the 
SU{N) Haldane- Shastry spin model with exchange interaction. It seems likely that 
our multi-component Hubbard model of the 1/r hopping is also completely integrable for 
arbitrary on-site energy U at arbitrary filling number |^ . However, we have not found any 



elegant way to obtain the wavefunctions and the energy spectrum for the finite on-site energy 
case. The nested Bethe-ansatz might be an essential key to the exact solutions at finite on- 
site energy. It is also very likely that the SU{N) system also exhibits a metal-insulator 
phase transition at half filling when changing the bandwidth and the on-site energy, as in 
the SU(2) case discovered by Gebhard and Ruckenstein about two years ago. It is also 
of great interest to study the ground state properties of the system as a function of the 
interaction strength, such as, the spin and charge susceptibilities and various ground state 
correlators. It also remains to find the integrability condition for the model at the finite 
on-site energy. 
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